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Abstract: This paper describes measurements of the sum of the transverse energy of 
particles as a function of particle pseudorapidity, rj, in proton-proton collisions at a centre- 
of-mass energy, y/s = 7 TeV using the ATLAS detector at the Large Hadron Collider. The 
measurements are performed in the region \rj\ < 4.8 for two event classes: those requiring 
the presence of particles with a low transverse momentum and those requiring particles 
with a significant transverse momentum. In the second dataset measurements are made in 
the region transverse to the hard scatter. The distributions are compared to the predictions 
of various Monte Carlo event generators, which generally tend to underestimate the amount 
of transverse energy at high \r)\. 
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1 Introduction 

The main aim of the Large Hadron Collider (LHC) general-purpose detectors is to explore 
physics in collisions around and above the electroweak symmetry-breaking scale. Such 
processes typically involve high momentum transfer, which distinguishes them from the 
dominant processes, namely low momentum transfer strong force interactions described 
by non-perturbative Quantum Chromodynamics (QCD). In order to collect enough data 
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to be sensitive to rare processes it is necessary to run the LHC at high instantaneous 
luminosities, meaning that multiple proton-proton interactions are very likely to occur 
each time the proton bunches collide. It is essential that the Monte Carlo event generators 
used to simulate these processes have an accurate description of the soft particle kinematics 
in inclusive proton-proton interactions over the entire acceptance of the LHC experiments, 
such that reliable comparisons can be made between theoretical predictions and the data 
for any process of interest. 

Protons are composite objects made up of partons, the longitudinal momentum dis- 
tributions of which are described by parton distribution functions (PDFs). When protons 
interact at the LHC the dominant parton-parton interaction is i-channel gluon exchange. 
Due to the composite nature of the protons it is possible that multiple parton-parton in- 
teractions (MPI) occur in the same proton-proton interaction. Therefore, if a hard parton- 
parton interaction occurs it will likely be accompanied by additional QCD interactions, 
again predominately low momentum t-channel gluon exchange. Any part of the interac- 
tion not attributed to the hard parton-parton scatter is collectively termed the underlying 
event, which includes MPI as well as soft particle production from the beam-beam rem- 
nants. Monte Carlo event generators that simulate any hard process at the LHC must also 
include an accurate description of the underlying event. 

At low momentum transfer, perturbative calculations in QCD are not meaningful and 
cross-sections cannot currently be computed from first principles. Phenomenological mod- 
els are therefore used to describe the kinematics of particle production in inclusive proton- 
proton interactions and in the underlying event in events with a hard scatter; these must 
be constrained by, and tuned to, data. 

This paper presents a measurement of the sum of the transverse energy, T,Et, of 
particles produced in proton-proton collisions at the LHC, using the ATLAS detector [1]. 
The SSx distribution is measured in bins of pseudorapidity 1 , rj, in the range < 4.8. 
Distributions of the S£"t and the function of \rj\ are presented. These 

measurements are performed with two distinct datasets. The first is as inclusive as possible, 
with minimal event selection applied, sufficient to ensure that an inelastic collision has 
occurred. This is termed the minimum bias dataset and is studied in order to probe the 
particle kinematics in inclusive proton-proton interactions. Understanding these processes 
is vital to ensure a good description of multiple proton-proton interactions in runs with high 
instantaneous luminosity. The second dataset requires the presence of two jets with high 
transverse energy, Et > 20 GeV, which ensures a hard parton-parton scatter has occurred 
and therefore allows the particle kinematics in the underlying event to be probed. This 
sample is termed the dijet dataset. Both datasets were collected during the first LHC runs 
at -y/s = 7 TeV in 2010. The data samples correspond to integrated luminosities of 7.1 fib^ 1 

1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the z-axis along the beam pipe. The i-axis points from the IP to the centre 
of the LHC ring, and the j/-axis points upward. Cylindrical coordinates (r, cj>) are used in the transverse 
(x — y) plane, <j> being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms 
of the polar angle 6 with respect to the beamline as r\ = — lntan(#/2). 
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for the minimum bias measurement 2 and 590 fib^ 1 for the dijet measurement. Such small 
data samples are used because the early LHC runs had a very low instantaneous luminosity 
ensuring a negligible contribution from multiple proton- proton interactions. The larger 
sample for the dijet analysis is used as the cross-section for a hard scatter is significantly 
lower than for inclusive proton-proton interactions. 

Many previous measurements of the kinematic properties of particles produced in 
minimum bias events [2-5] and in the underlying event [6-11] were restricted to the central 
region of the detectors. This is because they used tracking detectors, with limited coverage, 
to study charged particles, or because they used only the central region of the calorimeters, 
where the tracking detectors could be used for calibration. Measurements of the mean of 
the sum of the energy of particles as a function of in minimum bias events and in 
the underlying event were performed with the CMS forward calorimeter [12]; these were 
limited to the very forward region (3.15 < \r/\ < 4.9). LHCb has performed measurements 
of charged particle multiplicities in the regions —2.5 < r] < —2.0 and 2.0 < r] < 4.5 [13]. 

The measurements described in this paper utilize the entire acceptance of the ATLAS 
calorimeters, |?7| < 4.9, allowing the T,E^ to be probed and unfolded in the region \r]\ < 4.8. 
Unless otherwise stated, the central region will refer to the range \r]\ < 2.4 and the forward 
region will refer to the range 2.4 < \rj\ < 4.8. The measurement is performed with the 
ATLAS calorimeters and is corrected for detector effects so that the variables are defined 
at the particle- level (see Section 2), which includes all stable particles (those with a proper 
lifetime greater than 3 x 10 -11 seconds). Both the mean and distributions of the T,E^ are 
measured. This provides additional information, giving a complete picture of both inclusive 
proton-proton interactions and the underlying event in dijet processes, within the entire 
acceptance of the general purpose LHC detectors. The relative levels of particle production 
in the forward and central regions may be affected by the contribution from beam-beam 
remnant interactions, details of the hadronization as modelled with colour reconnection 
between quarks and gluons, the relative contribution from diffractive processes and the 
parton distribution functions in this kinematic domain. 

This paper is organized as follows. Section 2 defines the particle-level variables. Sec- 
tion 3 describes the Monte Carlo models that are used to correct the data for detector 
effects and to compare to the final unfolded results. The ATLAS detector is discussed in 
section 4, the event reconstruction in section 5 and the event selection in section 6. The 
method used to correct the data for detector effects is described in section 7. The sys- 
tematic uncertainties are described in section 8. Section 9 presents and discusses the final 
results and compares them to various Monte Carlo simulations. Finally, conclusions are 
given in section 10. 

2 Particle-level variable definitions 

In data, events are selected and variables defined using calibrated detector-level quantities. 
Corrections for detector effects are then applied. In order to compare the corrected data 
with predictions from Monte Carlo event generators without passing the events through 

2 The run dependence of the analysis was checked in a larger sample and found to be negligible. 
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a simulation of the ATLAS detector, it is necessary to define variables at the particle- 
level. The particle- level ££t is defined at the generator level by summing the E^ of all 
stable charged particles with momentum p > 500 MeV and all stable neutral particles with 
p > 200 MeV. Lower momentum particles are not included as they are unlikely to deposit 
significant energy in the ATLAS calorimeters. 

The S£/r distribution is defined as j^—^ , where N evt is the number of events in 
the sample. It is measured in six regions: 0.0 < \r)\ < 0.8, 0.8 < \rj\ < 1.6, 1.6 < \r]\ < 2.4, 
2.4 < \r/\ < 3.2, 3.2 < \r/\ < 4.0 and 4.0 < \r/\ < 4.8. In addition the mean T,E T over all 
events, per unit r]-cp, is measured as a function of \rj\. This is denoted as the transverse 
energy density (E^ nsity ) and is defined as ( ^?^ T )- I n the minimum bias measurement, 
the Y,Et includes particles at any cf). In the dijet measurement, the T,E^ is measured 
using only particles that are in the azimuthal region transverse to the hard scatter, namely 
^ < |A0| < ^p, where A(f> is the azimuthal separation between the leading jet and a given 
particle. This region of phase space contains limited particle production from the hard 
parton-parton interaction and is therefore most sensitive to the underlying event. 

2.1 Particle-level minimum bias event selection 

The events in the minimum bias analysis contain at least two charged particles with 
Pt > 250 MeV and \r/\ < 2.5, reflecting as closely as possible the requirement of a recon- 
structed vertex, as will be discussed in section 6. 

2.2 Particle-level dijet event selection 

The events in the dijet analysis contain at least two particle- level jets 3 . Both the leading 
and sub-leading jets must have £f > 20 GeV and \rf et \ < 2.5, reconstructed with the 
anti-kt [14] algorithm with radius parameter R = 0.4. This selection ensures that a hard 
scattering has occurred. A relatively small radius parameter reduces the probability of the 
jet algorithm collecting particles that are not associated with the hard scatter. In order to 
select a well balanced back-to-back dijet system, the jets satisfy |A(/>y| > 2.5 radians, where 

A0;i is the difference in azimuthal angle of the leading and sub-leading jet, and —J^t > 0.5, 

JJ E J T 

where is the E^ of the (sub-)leading jet. The latter requirement retains most of the 

dataset, but avoids topologies in which there is a large transverse energy difference between 
the leading and sub- leading jets. A well balanced dijet system suppresses contributions from 
multijet events, allowing a clearer distinction between regions with particle production 
dominated by the hard scatter and by the underlying event. 

3 Monte Carlo event generators 

This section describes the Monte Carlo event generator (MC) models used to correct the 
data for detector effects, to assign systematic uncertainties to the corrections due to the 
physics model, and for comparisons with the final unfolded data. The PYTHIA6 [15], 
PYTHIA8 [16], Herwig++ [17] and EPOS [18] generators are used, with various tunes that are 

3 A particle-level jet is built from all stable particles, excluding neutrinos and muons. 
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described below. First a brief introduction to the relevant parts of the event generators is 
given. 

PYTHIA 6 and PYTHIA 8 are general purpose generators that use the Lund string hadroniza- 
tion model [19]. In PYTHIA 6 there is an option to use a virtuality-ordered or pp-ordered 
parton shower, with the latter used in most recent tunes. In PYTHIA 8, the pp-ordered 
parton shower is used. The inclusive hadron-hadron interactions are described by a model 
that splits the total inelastic cross-section into non-diffractive processes, dominated by t- 
channel gluon exchange, and diffractive processes involving a colour-singlet exchange. The 
diffractive processes are further divided into single-diffractive dissociation, where one of 
the initial hadrons remains intact and the other is diffractively excited and dissociates, and 
double-diffractive dissociation where both hadrons dissociate. Such events tend to have 
large gaps in particle production at central rapidity. The smaller contribution from central 
diffraction, in which both hadrons remain intact and particles are produced in the central 
region, is neglected. The 2 — > 2 non-diffractive processes, including MPI, are described by 
lowest-order perturbative QCD with the divergence of the cross-section as pp — > regu- 
lated with a phenomenological model. There are many tunable parameters that control, 
among other things, the behaviour of this regularization, the matter distribution of par- 
tons within the hadrons, and colour reconnection. When pp-ordered parton showers are 
used, the MPI and parton shower are interleaved in one common sequence of decreasing 
pp values. For PYTHIA 6 the interleaving is between the initial-state shower and MPI only, 
while for PYTHIA 8 it also includes final-state showers. Since the pp-ordered showers and in- 
terleaving with MPI are considered to be a model improvement, the most recent PYTHIA 6 
tunes are made with this configuration. This is also the only configuration available in 
PYTHIA 8. A pomeron-based approach is used to describe diffractive events, using (by de- 
fault) the Schuler and Sjostrand [20] parameterization of the pomeron flux. In PYTHIA 6 
the diffractive dissociations are treated using the Lund string model, producing final-state 
particles with limited pp. In PYTHIA 8 the dissociations are treated like this only for events 
with a diffractive system with a very low mass; in higher mass systems diffractive parton 
distributions from HI [21] are used to include diffractive final states which are characteristic 
of hard partonic interactions. In this case, the full machinery of MPI and parton showers 
is used. This approach yields a significantly harder pp spectrum for final-state particles. 

Herwig++ is another general purpose generator, but with a different approach: it uses 
an angular-ordered parton shower and the cluster hadronization model [22]. It has an 
MPI model similar to the one used by the PYTHIA generators, with tunable parameters 
for regularizing the behaviour at very low momentum transfer, but does not include the 
interleaving with the parton showers. Inclusive hadron-hadron collisions are simulated 
by applying the MPI model to events with no hard scattering. It is therefore possible 
to generate an event with zero 2—7-2 partonic scatters, in which only beam remnants 
are produced, with nothing in between them. While Herwig++ has no explicit model 
for diffractive processes, these zero-scatter events will look similar to double-diffractive 
dissociation. 

EPOS is an event generator used primarily to simulate heavy ion and cosmic shower 
interactions, but which can also simulate proton-proton interactions. EPOS provides an 
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Generator 


Version 


lune 


PDh 


7 TeV data 










MB 


UE 


PYTHIA6 


6.423 


AMBT1 [24] 


MRSTLO* [25] 


yes 


no 


PYTHIA 6 


6.423 


DW [26] 


CTEQ5L [27] 


no 


no 


PYTHIA 6 


6.423 


PerugiaO [28] 


CTEQ 5L 


no 


no 


PYTHIA 8 


8.145 


4C [29] 


CTEQ6L1 [30] 


yes 


no 


Herwig++ 


2.5.1 


UE7-2 [31] 


MRSTLO** [25] 


no 


yes 



Table 1. MC tunes used to unfold the data and to determine the physics model dependent system- 
atic uncertainty. The last two columns indicate whether the data used in the tune included 7 TeV 
minimum bias (MB) and/or underlying event (UE) data. 



Generator 


Version 


Tune 


PDF 


7 TeV data 










MB UE 


PYTHIA 6 


6.425 


AUET2B : CTEQ6L1 [32] 


CTEQ 6L1 


no yes 


PYTHIA 8 


8.153 


A2:CTEQ6L1 [33] 


CTEQ 6L1 


yes no 


PYTHIA 8 


8.153 


A2:MSTW2008L0 [33] 


MSTW2008L0 [34] 


yes no 


EPOS 


1.99_v2965 


LHC 


N/A 


yes no 



Table 2. Additional MC tunes used to compare to the unfolded data only. The last two columns 
indicate whether the data used in the tune included 7 TeV minimum bias (MB) and/or underlying 
event (UE) data. 

implementation of a parton based Gribov-Regge [23] theory which is an effective, QCD- 
inspired field theory describing hard and soft scattering simultaneously. EPOS calculations 
thus do not rely on the standard PDFs as used in generators like PYTHIA and Herwig++. 
At high parton densities a hydrodynamic evolution of the initial state is calculated for the 
proton-proton scattering process as it would be for heavy ion interactions. The results 
presented here use the EPOS LHC tune, which contains a parameterized approximation of 
the hydrodynamic evolution. The optimal parameterization has been derived from tuning 
to LHC minimum bias data. 

The reference MC sample used throughout this study is the AMBT1 [24] tune of PYTHIA 6. 
In order to check the model dependence of the data corrections, additional generators and 
tunes are considered. These are summarized in table 1 along with information about the 
PDFs used and whether minimum bias or underlying event data at yfs = 7 TeV were used 
in the tune. Of the PYTHIA 6 tunes listed, only DW uses the old virtuality-ordered parton 
shower without interleaving with MPI. Some more recent tunes are also used to compare 
to the unfolded data; these are summarized in table 2. For these more recent tunes the 
PDF is explicitly given in the name as there are different instances of each tune that use 
different PDFs and hence have different parameters. 
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4 The ATLAS detector 



The ATLAS detector is described in detail in ref. [1]. Here only the components most 
relevant for this measurement are described. 

Tracks and interaction vertices are reconstructed with the inner detector tracking sys- 
tem, which consists of a silicon pixel detector, a silicon strip detector and a transition ra- 
diation tracker, all immersed in a 2 T axial magnetic field. The calorimeter systems are of 
particular importance for the measurements presented in this paper. The ATLAS calorime- 
ter system provides fine-grained measurements of shower energy depositions over a wide 
range of ij. A highly segmented electromagnetic liquid argon (LAr) sampling calorimeter 
covers the region < 3.2, with granularity that ranges from 0.003x0.10 or 0.025x0.025 
to 0.1x0.1 in At] x A<j>, depending on depth segment and pseudorapidity. It is divided 
into a barrel part (\r}\ < 1.475) and an endcap part (1.375 < \rj\ < 3.2). The hadronic 
barrel (\rj\ < 1.7) calorimeter consists of steel absorbers and active scintillating tiles, with 
a granularity of either 0.1x0.1 or 0.2x0.1 depending on the layer. The hadronic endcap 
(1.5 < < 3.2) and forward (3.1 < \r}\ < 4.9) electromagnetic and hadronic calorimeters 
use liquid argon technology. The granularity in the hadronic endcap ranges from 0.1x0.1 
to 0.2x0.2. In the forward calorimeter, the cells are not arranged in projective towers but 
are aligned parallel to the beam axis. As such the readout granularity is not constant in 

v~4>- 

Minimum bias trigger scintillator (MBTS) detectors are mounted in front of the endcap 
calorimeters on both sides of the interaction point and cover the region 2.1 < < 3.8. The 
MBTS is divided into inner and outer rings, both of which have eight-fold segmentation, 
and is used to trigger the events analysed in this paper. 

5 Event reconstruction 

This analysis is based on topological clusters in the calorimeter, which represent an at- 
tempt to reconstruct three-dimensional energy depositions associated with individual par- 
ticles [35]. The topological clustering algorithm proceeds through the following steps. First, 
seed cells are found that have \E\ > 4a above the noise level, where E is the cell energy 
measured at the electromagnetic scale 4 and calibrated using test-beam data [36-39]. Next, 
neighbouring cells are collected into the cluster if they have \E\ > 2a above the noise level. 
Finally, all surrounding cells are added to the cluster until no further cells with \E\ > 2a 
are among the direct neighbours. 

The detector-level T,Et is formed by summing the Et of all clusters in the r]-(j) region 
of interest. Negative energy clusters are included, leading to a convenient cancellation of 
the contributions from noise, which can be either negative or positive. 

To correct these clusters back to the particle-level, it is first necessary to determine the 
particle momenta to which the ATLAS calorimeters are sensitive. Using a GEANT4 [40] sim- 

4 The electromagnetic scale is the basic calorimeter signal scale for the ATLAS calorimeters. It gives the 
correct response for the energy deposited in electromagnetic showers, but does not account for the lower 
response to hadrons. 
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Figure 1. The di-photon invariant mass in the region (a) 1.52 < 77 < 2.37 and (b) 4.2 < 77 < 4.8. 
The data are compared to the MC simulation with the best fit scale factor applied (this is 0.97±0.02 
for (a) and 1.01±0.02 for (b)). The contribution from the MC signal ir° — > 77 templates and 
background templates are also shown separately. The arrows indicate the fit range. 

ulation of the ATLAS detector [41], generator-level particles are propagated from the pri- 
mary vertex to the calorimeters and the fraction of their energy deposited in the calorime- 
ters as clusters is studied as a function of \rj\ and p of the particle. As discussed in section 2, 
charged particles with p > 500 MeV and neutral particles with p > 200 MeV are found to 
deposit enough energy in the calorimeter to be included in the particle-level definition for 
all \i]\ regions. Particles with lower momenta contribute a negligible amount to the cluster 
X_E/r and are therefore excluded from the particle-level T,Ej< definition. 

In order to properly correct for detector effects, the detector simulation must accurately 
describe the energy response of the calorimeters to low energy particles. The simulation 
calibration is refined using the di-photon invariant mass distribution of 7r° — > 77 candidates. 
In data selected with the MBTS trigger, pairs of photon candidates in a given 77 region are 
formed and their invariant mass, m 77 , is constructed. In order to reduce combinatorial 
background, only events with exactly one pair in the 77 region are considered. The data are 
compared to MC signal plus background templates in 77 bins, which are chosen to reflect 
the boundaries of the calorimeter sub-systems. 

The signal templates are derived from the PYTHIA 6 AMBT1 samples, by matching pairs 
of clusters to generator level photons from a 7r° decay. The background templates are 
obtained using pairs of clusters that are not matched. The energies of the clusters in the 
signal template are scaled by an energy response scale factor. This is varied and the \ 2 
between the data and MC distributions is minimized in order to determine the best fit 
value. Deviations from unity are typically 2-3% but reach values of up to 10% in some 77 
regions. This scale factor is then applied to the energy of the MC clusters before unfolding 
the data. Figure 1 shows the m 77 distribution in data compared to the MC in two sample 
1 77 1 regions with the best fit scale factor applied. 
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6 Event selection 



Events in the minimum bias analysis are selected with a one-sided MBTS trigger, which 
requires one counter on either side of the detector to be above noise threshold, suppressing 
contributions from empty beam crossings and beam-induced background. In order to 
suppress these contributions further, events are required to have a reconstructed primary 
vertex with at least two associated tracks with pp > 150 MeV and |7y| < 2.5. Note that the 
track px cut is lower than the 250 MeV particle-level cut described in section 2.1. This 
is because tracking and vertex reconstruction inefficiencies result in events with at least 
two 150 MeV reconstructed tracks having the same £'^ ensity as events with at least two 
250 MeV charged particles, according to the MC models considered in this analysis. 

Furthermore, events having more than one reconstructed vertex with five or more 
tracks are vetoed to suppress contributions from multiple proton-proton interactions. Five 
tracks are required on the additional vertices so that events with secondary vertices from 
decaying particles are not vetoed. 

Events in the dijet analysis are also selected with the one-sided MBTS trigger and 
are required to pass the same event selection criteria as the minimum bias analysis. In 
addition, they are required to contain two back-to-back jets passing the same kinematic 
selection criteria as the particle-level jets described in section 2.2. 

7 Corrections for detector effects 

The SE't distributions are unfolded in each \rj\ region using an iterative Bayesian unfolding 
technique [42]. The E^ ensity distribution is obtained by taking the mean of each unfolded 
distribution and dividing by the \r}\ and <t> phase space. An unfolding matrix is 
formed from events generated with PYTHIA 6 AMBT1, passed through the GEANT4 simulation 
of the ATLAS detector. The detector simulation accounts for energy losses of the particles 
in material upstream of the calorimeter, for charged particles that bend in the magnetic 
field and get swept out of the calorimeter acceptance, and for the calorimeter response 
and resolution. Before unfolding each EEr distribution, the MC is reweighted by a fit 
to the ratio of the data to the MC detector-level SEp distribution, so that the T,E^ 
distribution matches that seen in data. The MC significantly underestimates the T,E^ in 
the forward region, as seen in figure 2, where the detector-level T,E^ distribution in the 
region 4.0 < \rj\ < 4.8 is shown for both data and MC, before and after reweighting, for 
both the minimum bias and dijet selections. 

The unfolding matrix associates the T,E^ formed from clusters with the EEr formed 
from generator-level particles. Events that pass the detector-level but not the particle-level 
selection criteria and vice versa are also accounted for in the correction procedure. The prior 
distribution of the particle-level EEt is initially taken from PYTHIA 6 AMBT1 (reweighted 
to data) and the unfolding procedure is iterated twice, with the prior distribution replaced 
by the unfolded distribution after each iteration. A stable result is achieved after two 
iterations. 
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Figure 2. The detector-level ££>r distribution in the region 4.0 < \r]\ < 4.8 for data compared to 
the nominal detect or- level PYTHIA 6 AMBT1 prediction and the reweighted detect or- level PYTHIA 6 
AMBT1 prediction in (a) the minimum bias events and (b) the dijet events. 

8 Systematic uncertainties 

The dominant systematic uncertainties arise from three sources: (1) the accuracy with 
which the MC simulates the energy response of the calorimeters to low energy particles, 
(2) the knowledge of the amount of material upstream of the calorimeters and (3) the MC 
generator model dependence in the unfolding. In the dijet analysis an additional uncer- 
tainty arises from the accuracy with which the MC simulates the jet energy scale. These 
sources are discussed in the following sub-sections. In each case the uncertainty on the 
unfolded data is obtained by shifting the MC by ±1ct for the source in question and com- 
paring with the nominal unfolded data. In order to give information about the correlations 
of the systematic uncertainties between bins and between the different distributions in this 
paper, each source is split into different components. These systematic uncertainties are 
summarized in tabular form in appendix A. 

The following additional potential sources of systematic uncertainty are found to be 
negligible: energy resolution, multiple proton-proton interactions, contributions from noise 
and beam-induced backgrounds, simulation of the primary vertex position, simulation of 
the trigger selection, and simulation of the position of the forward calorimeter. 

8.1 Calorimeter energy response 

The systematic uncertainty on the calorimeter energy response is determined separately for 
electromagnetic and hadronic particles. An average is then obtained, using the PYTHIA 6 
AMBT1 prediction of the relative contributions to the by different particle types. For 
electromagnetic particles the systematic error comes from uncertainties on the extraction 
of the energy scale from fits to the m 77 distributions in 7r° — > 77 candidates. These are 
obtained from variations in the fit range, the background shape, the criteria for matching 
reconstructed photons to generator-level photons in the production of the signal template, 
variations in the simulation of the calorimeter resolution, and consistency with a similar 
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analysis using tighter kinematic and photon identification cuts. The total uncertainty 
depends on the \rj\ region and is generally at the level of 2-4%, but increases up to 15% in 
the regions where different calorimeter sub-systems overlap. 

The uncertainty on the energy response for hadronic particles in the central region, 
where there is good coverage from the inner tracking detector, is obtained from studies of 
the ratio of the calorimeter energy measurement to the inner detector track momentum 
measurement, for isolated charged pions [43]. The uncertainty is obtained by taking the 
difference between data and MC in p and \r]\ bins and is found to be 3.5% for \r/\ < 0.8 
and 5% for 0.8 < < 2.4. In the forward region the energy response uncertainty for 
hadrons is taken from the difference between the MC and data in test-beam studies of 
charged pions [44]. This leads to a one-sided uncertainty for hadrons relative to electro- 
magnetic particles of +5% in the region 2.5 < \rj\ < 3.2 and +9% in the forward calorimeter 
(\V\ > 3.2). 

The only component of the systematic uncertainty on the energy response assumed to 
be correlated between \r]\ bins is that of the forward calorimeter (determined from test- 
beam results), which affects the bins 3.2 < |r/| < 4.0 and 4.0 < \r/\ < 4.8. The systematic 
uncertainties from the 7r° — > 77 fits are assumed to be uncorrelated as the m 77 shapes are 
rather different in the different \rj\ regions, resulting in different possible systematic shifts. 
Similarly, the difference between data and MC for the ratio of calorimeter energy to inner 
detector track momentum does not show systematic shifts in one direction and is assumed 
to be uncorrelated. 

Appendix A gives both the uncorrelated and correlated uncertainties in each bin of 
each distribution. The former vary between 2.4% and 5.4% for the E^ ensity in the min- 
imum bias data, depending on the |?7| region. The largest uncertainty is in the region 
0.8 < \rj\ < 1.6, which contains the region of overlap between the barrel and endcap elec- 
tromagnetic calorimeters (1.375 < [77! < 1.475). The correlated source is about —6% for 
the two highest \rj\ bins in the minimum bias data and about —8% in the dijet data. Note 
that a positive uncertainty on the energy scale in the MC leads to a negative uncertainty 
on the corrected result in the data. The uncertainty is higher in the dijet data, due to a 
larger contribution from events where the detector-level jets pass the selection criteria but 
the generator- level jets do not. Their Y,E<i distribution is taken from the MC so a shift in 
the energy scale leads to an additional bias in the corrected result. 

8.2 Material description 

The amount of material upstream of the calorimeters affects the T,E^ distributions because 
particles can interact and lose some of their energy before reaching the calorimeter. It is 
therefore important to have a realistic description of the material in the MC simulation 
used to perform the detector corrections. 

In order to assess the systematic uncertainty arising from possible discrepancies in 
the material description, detector corrections are recalculated using a special PYTHIA6 
AMBT1 sample with additional material. The sample is based on a similar one described 
in section 3 of ref. [45], but with additional material introduced in the forward region. 
The results are compared to the nominal unfolded data and the difference is taken as a 
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symmetric systematic uncertainty to account for the possibility of the MC simulation either 
underestimating or overestimating the amount of material. 

In order to understand the correlations between the uncertainties in different bins, 
the additional material is split into three components: (1) extra material upstream of 
the barrel calorimeter, (2) an increase in material density in the barrel-endcap overlap 
region and (3) additional material in the inner detector, the inner detector services and the 
forward region, as well as an increase in the material density in some detector volumes in 
the forward region. The systematic uncertainties arising from components (1) and (2) are 
assumed to be correlated between \r]\ bins, whereas the uncertainty arising from component 
(3) is assumed to be uncorrelated, due to the fine structure of these detectors with respect 
to the wide bins used in this analysis. 

Source (1) affects only the first two \rj\ bins, at the level of about 3% in the minimum 
bias data and 1.3-2.5% in the dijet data. The uncertainty is generally smaller in the dijet 
data as the particles in these events tend to have larger momenta. Source (2) affects only 
the second and third bins and is less than 1%. Source (3) affects all \r]\ bins and ranges 
between 0.23% and 5.5%, with the largest uncertainty in the region 1.6 < \r]\ < 2.4, where 
there is a large amount of material associated with the inner detector. 

8.3 Physics model dependence 

The MC model used to correct the data can affect the results as a realistic description 
of particle kinematics is needed. The model dependence is minimized by first reweighting 
the detector-level MC to the data and then by iterating the unfolding, using the unfolded 
data as the new prior distribution after each iteration. This reduces the dependence on 
the EE'x spectrum itself; however, other kinematic distributions can also affect the un- 
folding. One important variable is the -E/p of the individual particles, as the calorimeter 
response to a particle is energy dependent. The dependence on the model is investigated 
by performing the unfolding with other MC models. The following MC models and tunes 
are considered: PYTHIA6 AMBT1 (nominal), PYTHIA6 DW, PYTHIA6 PerugiaO, PYTHIA8 4C 
and Herwig++ UE7-2. Details of these tunes are given in table 1. The MC model used to 
assess the systematic uncertainty is chosen to ensure a reasonable spread in the particle 
kinematics with respect to the reference PYTHIA 6 AMBT1 model. Figure 3 shows distribu- 
tions of -^t3t x d\E T \ ' wnere -E"T 0t * s the sum over even ts of the detector- level E-E/r> and 
-E/p is the detector-level cluster transverse energy. These distributions show the relative 
contribution to the TiE^ from clusters with a given E^£. \E^\ is plotted instead of E^ since 
the former leads to a cancellation in the contribution from noise. Figure 3(a) shows the 
distribution in minimum bias events for the region 3.2 < < 4.0. This region is shown 
as it has significant differences between data and MC. The contribution to the E_Et from 
high E^y clusters is smaller in data than in PYTHIA 6 AMBT1. The model with the largest 
deviations from PYTHIA 6 AMBT1 is Herwig++ UE7-2, indicating that this model can be used 
to assess possible biases in the unfolding due to this effect. It should be noted that at 
high |£t| Herwig++ UE7-2 lies above PYTHIA 6 AMBT1 while the data lie below it, but the 
final systematic uncertainty is symmetrized. The same distribution is shown in figure 3(b) 
for the sub-sample of events with £.&r > 15 GeV. Again, the data have a softer cluster 
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Figure 3. Distribution of x d\E T \ > wnere ^t * ^ s * ne sum over even ts of the detector-level 
Si?T, and Et is the detector-level cluster transverse energy (a) in minimum bias events in the 
region 3.2 < \rj\ < 4.0; (b) as in (a) but for events with T,Et > 15 GeV; (c) in dijet events in the 
region 0.0 < |?j| < 0.8; (d) as in (c) but for events with T,Et > 15 GeV. The data are compared to 
various MC predictions. 

|.&r| distribution. Here PYTHIA6 DW shows the largest deviations from PYTHIA6 AMBT1. 
Since unfolding with PYTHIA6 DW results in a larger shift in the corrected data than un- 
folding with Herwig++ UE7-2, the former is used to assess the systematic uncertainty in 
the minimum bias events. 

Figure 3(c) shows the same distribution in dijet events; the most central region is 
shown as the differences between data and MC are largest in this region. This time the 
data distribution is harder than PYTHIA6 AMBT1. Again Herwig++ UE7-2 has the largest 
deviations. Figure 3(d) shows the same distribution for events with Y,Et > 15 GeV; all 
the models agree well with the data, but Herwig++ UE7-2 has the largest deviations. For 
the dijet selection Herwig++ UE7-2 is therefore used to assess the systematic uncertainty. 

For both the minimum bias and dijet analyses, this systematic uncertainty is sym- 
metrized and treated as correlated between \rj\ bins (although not correlated between the 
two analyses). The uncertainties on the £^ ensity range from 2-4% for the minimum bias 
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data and are 2% or less for the dijet data. 
8.4 Jet energy scale 

In the dijet selection, events are required to contain at least two jets with > 20 GeV. 
It is possible that events that satisfy the detector-level criteria do not satisfy the criteria 
at the particle-level, and vice versa. This is accounted for in the correction procedure, 
but if there are differences in the jet energy scale between data and MC simulation this 
could result in a bias in the correction procedure. The uncertainty on the jet energy scale 
is described in ref. [46]. The corresponding uncertainty on the E^ nsity is at the level of 
1.6% in the most central bin and decreases to 0.13% in the most forward bin. It is treated 
as correlated between |r/| bins. For the HE^ distributions this source of uncertainty is 
negligible and therefore neglected in the region \rj\ > 2.4. 



9 Results 

9.1 Nominal Results 

The unfolded E^ ensity distributions are shown in figure 4 for both the minimum bias and 
the dijet selections. The filled bands indicate the systematic and statistical uncertainties 
on the data, added in quadrature. In all bins the systematic uncertainty is significantly 
larger than the statistical uncertainty. The E^ ensity distribution in the minimum bias data 
dips in the central region. Since the relative fraction of low momentum particles is higher 
in the central region than in the forward region, fewer central particles pass the selection 
criteria described in section 2, hence reducing the T^Et m the central region. The dip in 
the central region is less prominent in the dijet data; this feature is discussed below. 

Figure 5 shows the ratio of the E^ ensity in the dijet transverse region to the £'^ ensity in 
minimum bias events. The correlations between the systematic uncertainties for the dijet 
and minimum bias distributions are taken into account. All systematic uncertainties but 
the physics model dependence and jet energy scale are taken as correlated between the 
two. The £^ enhlty in the transverse region for the dijet selection is larger than the E^ ensity 
in the minimum bias data. This increase is expected, due primarily to the presence of a 
hard scatter, which will bias the selected events away from peripheral proton scatters and 
towards small impact parameter ("head-on") proton- proton interactions. This means that 
more parton-parton interactions are likely to occur in the underlying event in the dijet 
data than in the collisions with a larger impact parameter that characterize the events in 
the minimum bias dataset. 

The unfolded data are compared to various MC models. In the minimum bias sample 
the E^ ensity distribution in figure 4(a) is well described by PYTHIA 6 AMBT1 in the central 
region. This is expected as this tune was prepared with ATLAS 7 TeV minimum bias data in 
the region \r]\ <2.5 [3]. At higher \rj\ values, however, the E^ nsity is underestimated and is 
approximately 25% too low in the highest \r]\ bin. The PYTHIA 6 AUET2B : CTEQ6L1 prediction 
is very similar to that from PYTHIA 6 AMBT1, with slightly more energy in the central region 
and less in the forward region, meaning that the description of the |ry| dependence is even 
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Figure 4. Unfolded E T cnsl y distribution compared to various MC models and tunes for (a) the 
minimum bias selection and (b) the dijet selection in the transverse region. The filled band rep- 
resents the total uncertainty on the unfolded data. N c h refers to the number of charged particles 
in the event, and and r] ch are, respectively, the pt and 77 of those particles. 2Vj B t refers to the 
number of jets, E^ 11 ^ is the Ei of the (sub-)leading jet, rf ct is the jet pseudorapidity, and A^jj is 
the azimuthal angle difference between the two leading jets. p ch ( ncutral ) refers to the momentum of 
the charged(neutral) particles used in the T,E^ calculation. 
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Figure 5. Unfolded E^ enslty distribution in the dijet data transverse region divided by that in the 
minimum bias data, compared to various MC models and tunes. The filled band represents the 
total uncertainty on the unfolded data. 



worse. PYTHIA6 DW underestimates the E^ ensity in all \rj\ bins. Despite this it provides 
an improved description of the \rj\ dependence of the £'^ ensity . PYTHIA8 4C overestimates 
the E^ ensity in the central region. The agreement improves in the region 1.6 < |r/| < 3.2, 
but in the higher |r/| bins the £'^ ensity is underestimated. Herwig++ UE7-2 overestimates 
the E^ nsity in the central region, describes the data well in the region 2.4 < \r/\ < 3.2, and 
undershoots the data at higher \rj\. The EPOS LHC prediction provides the best description 
over the entire \rj\ region, although it does fall slightly too fast with \r]\. It should be 
noted that, with the exception of EPOS LHC and PYTHIA 6 DW, while some models and tunes 
appear to agree better in some regions than others, this is generally due to differences in 
the total level of particle production. The overall pattern remains the same: the ^ ensity 
in the forward region is too low relative to the central region. 

In the dijet selection in figure 4(b), all of the MC models and tunes perform reasonably 
well in the central region, apart from EPOS LHC which underestimates the E^? nsity in all \rj\ 
bins. PYTHIA 6 AUET2B : CTEQ6L1 slightly overestimates the energy in the most central bins, 
and all the other predictions are slightly too low. As was the case in the minimum bias 
analysis, the £^ ensity in the forward region is underestimated. PYTHIA 8 4C is approximately 
20% too low in the most forward bin, while PYTHIA 6 AMBT1, Herwig++ UE7-2 and PYTHIA 6 
AUET2B : CTEQ6L1 are 25-30% too low. PYTHIA 6 DW provides the best description of the \rj\ 
dependence, although the overall amount of energy is too low. 

The fall-off with \rj\ of the ratio of the E^ ensity in dijet and minimum bias events seen in 
figure 5 is reproduced by the models, with PYTHIA 6 AMBT1 and AUET2B : CTEQ6L1 describing 
the data the best. The reduction in the ratio with \r/\ is partly due to the momentum cuts 
on the particles included in the SSt calculation. In the dijet data, the particles tend to 
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have larger momenta and so fewer are removed from the calculation. According to 
PYTHIA6 AMBT1, the momentum cuts remove 25(18)% of the .E^, ensity in the most central 
bin and a negligble amount in the most forward bin for the minimum bias (dijet) selections. 
The PYTHIA6 AMBT1 (no p cuts) curve in figure 5, shows the ratio when the momentum 
cuts on the particles contributing to the SE't have been removed. There is still a residual 
decrease with \rj\ which may be due to a contribution to the underlying event in the central 
region coming from particles associated with the hard scatter. 

The unfolded T,Eic distributions are shown in figures 6 and 7 for the minimum bias 
and dijet selections, respectively. The distribution peaks at higher values of E£t in the 
forward region due to the particle momentum cuts discussed above. In the region \rj\ < 3.2 
the distribution is broader than in the forward region, with more events populating the 
high X£t tails. There is therefore more event-by-event variation in the £2?t in the central 
part of the detector. These features are reproduced by the MC predictions. PYTHIA 6 AMBT1 
provides the best description of the ££t shape in the central region for the minimum bias 
data. For the dijet data, most of the tunes do a reasonable job, although PYTHIA 8 4C and 
EPOS LHC underestimate the high SE'x tails. As with the £^ ensity distributions, the X&r 
in the forward region is underestimated for all but the dijet PYTHIA 6 DW prediction. 

In summary, all of the MCs underestimate the amount of energy in the forward region 
relative to the central region, in both the minimum bias data and the underlying event, 
with the exception of PYTHIA 6 DW which provides a reasonable description of the dijet data, 
although the prediction is approximately one standard deviation below the central values 
measured in the data in all \rj\ bins. EPOS LHC provides the best overall description of 
the minimum bias data. PYTHIA 6 AMBT1 provides the best description in the most central 
region (\rj \ < 1.6), while at higher \r]\ values PYTHIA 8 4C and Herwig++ UE7-2 reflect the 
data more accurately. In the dijet analysis, all the MCs provide a reasonable description 
in the central region, apart from EPOS LHC. 

9.2 Variation in diffractive contributions 

In order to investigate the sensitivity of the £^ ensity to the fraction of diffractive events, fig- 
ure 8 compares the unfolded E^ ensity distribution in the minimum bias data to PYTHIA 8 4C 
with the nominal diffractive cross-sections (50.9 mb, 12.4 mb and 8.1 mb for non-diffractive, 
single-diffractive and double-diffractive processes, respectively) and to samples where the 
diffractive cross-sections have been doubled or halved, with the non-diffractive cross-section 
held constant. This is achieved by combining the separate MC samples for the different 
processes with adjusted weights, rather than by changing the relevant parameters when 
generating the samples. Diffractive processes tend to have less particle production than 
non-diffractive processes. As expected, increasing the diffractive contribution decreases the 
-^density. However, the shape of the _E^, ensity distribution is not significantly affected. 

9.3 Variation in parton distribution functions 

The overall energy as well as its \rj\ dependence are affected by the PDFs used as input to 
the MC model. In order to investigate the dependence on the PDFs, comparisons are made 
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Figure 6. Unfolded T.E^ distributions compared to various MC models and tunes for the minimum 
bias selection in the following regions: (a) 0.0 < \r)\ < 0.8, (b) 0.8 < |ry| < 1.6, (c) 1.6 < \t}\ < 2.4, 
(d) 2.4 < \r)\ < 3.2, (e) 3.2 < \rj\ < 4.0 and (f) 4.0 < \r}\ < 4.8. The filled band in each plot represents 
the total uncertainty on the unfolded data. A^h refers to the number of charged particles in the 
event, and and r/ ch are, respectively, the px and 77 of those particles. p ch ( noutra i) refers to the 
momentum of the charged (neutral) particles used in the Ei?x calculation. 
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Figure 7. Unfolded YiE^ distributions compared to various MC models and tunes for the dijet 
selection in the transverse region in the following \rj\ regions: (a) 0.0 < \r/\ < 0.8, (b) 0.8 < \r)\ < 1.6, 
(c) 1.6 < \r)\ < 2.4, (d) 2.4 < \t]\ < 3.2, (e) 3.2 < \r)\ < 4.0 and (f) 4.0 < \r]\ < 4.8. The filled band 
in each plot represents the total uncertainty on the unfolded data. Nj ct refers to the number of 
jets, E^ 11 ^ is the Et of the (sub-)leading jet, 7^°* is the jet pseudorapidity, and A0jj is the 
azimuthal angle difference between the two leading jets. p ch ( neutral ) refers to the momentum of the 
charged(neutral) particles used in the E£/r calculation. 
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Figure 8. Final unfolded i5 T ° nM y distribution for the minimum bias selection compared to PYTHIA 8 
4C with the nominal diffractive cross-sections, as well as enhanced and suppressed diffractive cross- 
sections, as described in the text. The filled band represents the total uncertainty on the unfolded 
data. Nch refers to the number of charged particles in the event, and p^ and r] ch are, respectively, 
the pt and r] of those particles. p ch ( neutral ) refers to the momentum of the charged(neutral) particles 
used in the T,E T calculation. 

between the data and the PYTHIA 8 A2 family of tunes, which use different input PDFs [33], 
with the following variations: 

1. Tune A2:CTEQ6L1. 

2. The A2:CTEQ6L1 tune parameters, but with the MSTW2008L0 PDFs. 

3. Tune A2:MSTW2008LD. 

4. Tune A2:CTEQ6L1 where the E^ ensity has been scaled by 0.93(0.96) for the minimum 
bias (dijet) selection so that it matches A2 :MSTW2008L0 in the most central bin. 

These comparisons are shown in figure 9. The first thing to note is that moving from 
the CTEQ6L1 to the MSTW2008L0 PDFs (and keeping all tune parameters the same) de- 
creases the amount of energy in the central region, but increases it in the forward region, 
presumably due to the increase in both the high-x and low- a; gluon PDF with respect to 
the mid-x region, where x is the proton momentum fraction carried by the gluon. When 
the parameters are tuned to data in the central region, the energy increases for the mini- 
mum bias prediction. If the £^ ensity obtained using A2:CTEQ6L1 is scaled down to match 
A2 : MSTW2008L0 in the most central bin, it is clear that the latter provides a better descrip- 
tion of the data in the forward region, with the underestimation in the most forward bin 
improving from about 30% to 15%. 
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Figure 9. Final unfolded _E T ° nsl y distribution compared to PYTHIA8 with variations of the PDFs 
used, as discussed in the text for (a) the minimum bias selection and (b) the dijet selection. The filled 
band represents the total uncertainty on the unfolded data. JV c h refers to the number of charged 
particles in the event, and p^jJ 1 and r] ch are, respectively, the px and r\ of those particles. iVj e t refers 
to the number of jets, E^ 11 ^ is the Ei of the (sub-)leading jet, rf ct is the jet pseudorapidity, 
and A</>jj is the azimuthal angle difference between the two leading jets. p ch ( neutral ) refers to the 
momentum of the charged (neutral) particles used in the ££/r calculation. 
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10 Conclusions 



Measurements of the E T ensi y and the SEt distributions as functions of \rj\ have been 
presented for two event classes: those requiring the presence of particles with a low trans- 
verse momentum (minimum bias) and those requiring particles with a significant transverse 
momentum (dijets), using proton-proton collision data at y/s =7 TeV recorded by the AT- 
LAS detector. In the dijet selection the distributions are measured in the region transverse 
in (f) to the hard scatter, in order to probe the particle production from the underlying 
event. The measurements are performed in the region |r/| < 4.8 for charged particles with 
p > 500 MeV and neutral particles with p > 200 MeV, and are the first to utilize the en- 
tire acceptance of the ATLAS calorimeters to probe the overall properties of inclusive 
proton-proton collisions, as well as the underlying event. The distributions are compared 
to various MC models and tunes. In general all MC predictions are found to underestimate 
the amount of energy in the forward region relative to the central region by 20-30%, with 
the exception of the PYTHIA 6 DW tune and EPOS LHC for the minimum bias data, although 
PYTHIA6 DW underpredicts the overall energy by 20-30%. For the PYTHIA 8 A2 tune series, 
this is improved if the MSTW2008L0 PDFs are used instead of the CTEQ 6L1 PDFs. 
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A Tabulated results and uncertainties 



The unfolded data are presented in tabular form in this appendix for the 2? T ensi y and the 
six SSt distributions, for both the minimum bias and dijet selections. Tables 3, 4 and 5 
give the unfolded data and systematic uncertainties for the £^ ensity for the minimum bias 
selection, the dijet selection, and the ratio between them, respectively. Tables 6-11 give the 
unfolded data and systematic uncertainties for the T,E^ distributions for the minimum bias 
selection and tables 12-17 give the corresponding information for the dijet analysis. In each 
case, the breakdown of the systematic uncertainties by source is also given. Each systematic 
source is described in section 8. The uncorrelated calorimeter energy scale systematic is 
denoted as E^' b > c ' d ' e 'f f or each of the six \rj\ regions, respectively. The correlated calorimeter 
energy scale systematic is denoted as Ei- The two correlated material systematic sources 
are denoted as M\ and M2 and the uncorrelated source is denoted as M^ ,b,c,d,e ^ for the six 
I rj I regions. All the above sources are correlated between the minimum bias data and the 
dijet data and therefore have the same symbol. 

The physics model systematic uncertainty on the minimum bias and dijet results, and 
on their ratio, are denoted as, P\, P2 and P, respectively. The jet energy scale systematic 
uncertainty is denoted as J. The physics model and jet energy scale systematic sources 
are uncorrelated between the minimum bias and dijet data. For the distributions J 
is negligible and therefore neglected in the region \rj\ > 2.4. 

The correlations between bins of a given distribution are indicated by the sign of the 
uncertainty. For example, in table 6 the uncertainty Ef is ± in the first three bins and =p 
in the remaining bins. This means that the first three bins are correlated with each other 
and anti-correlated with the remaining bins (a downward shift in the 'EEt will shift the 
low ££t bins up and the high T,E^ bins down). Since the individual sources within a 
given distribution are uncorrelated, the relationship between ± and =F between sources is 
not relevant to the calculation of the total error in a given bin. 

The uncertainties are given to two significant figures or a precision of 0.01%, whichever 
is smaller. In cases where the + and — uncertainty have a different precision the lowest 
precision is chosen for both. In cases where the uncertainty is not applicable, this is 
indicated with a dash. 
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\v\ 




\ dnrirh 1 


Stat. 


E{ 


E2 


Mi 


M 2 


M 3 * 


Pi 




Total 






[GeV] 


[%] 


[%} 


[%] 


[%} 


[%] 


[%] 


[%} 




[%] 


0.0 


- 0.8 


0.753 


±0.19 


+3.2 
-2.9 




±2.9 




±0.51 


±2. 


.6 


+5.1 
-4.9 


0.8 


- 1.6 


0.844 


±0.17 


+5.4 
-4.9 




±3.2 


±0.49 


±1.2 


±4. 


.6 


+7.9 
-7.5 


1.6 


- 2.4 


0.902 


±0.16 


+4.0 
-3.8 






±0.89 


±5.0 


±3. 


.4 


+7.4 
-7.2 


2.4 


- 3.2 


0.932 


±0.16 


+2.4 
-5.0 








±3.0 


±2. 


.5 


+4.6 

-6.4 


3.2 


- 4.0 


0.850 


±0.15 


+4.3 
-4.4 


-6.2 






±2.7 


±3. 


.2 


+6.0 
-8.7 


4.0 


- 4.8 


0.750 


±0.14 


+2.7 
-2.7 


-6.8 






±0.8 


±3. 


.6 


+4.6 
-8.2 



Table 3. Measured E T cnsl y and systematic uncertainty breakdown for the minimum bias data. 
The systematic uncertainties marked with a * are uncorrelated between \r)\ bins. 



\v\ 




id z T,E T \ 
\ drjd<f> ' 


Stat. 


E\ 


E 2 Mi 


M 2 


M 3 * 


P2 




J 


Total 






[GeV] 


[%] 


[%} 


[%] [%] 


[%] 


[%] 


[%] 




[%} 


[%] 


0.0 


- 0.8 


2.22 


±0.61 


+4.3 
-4.2 


±1.3 




±0.23 


±2. 


2 


+1.6 
-1.3 


+5.3 
-5.1 


0.8 


- 1.6 


2.37 


±0.54 


+7.2 
-6.4 


±2.5 


±0.38 


±0.96 


±0. 


12 


+1.3 
-1.3 


+7.8 
-7.1 


1.6 


- 2.4 


2.35 


±0.52 


+5.3 
-5.0 




±0.97 


±5.5 


±0. 


,41 


+0.98 
-0.92 


+7.8 
-7.6 


2.4 


- 3.2 


2.27 


±0.50 


+3.8 
-7.0 






±0.64 


±0. 


.55 


+0.80 
-0.37 


+4.0 
-7.1 


3.2 


-4.0 


1.88 


±0.51 


+6.1 
-5.8 


-8.2 




±1.1 


±1. 


3 


+0.46 
-0.17 


+6 
-10 


4.0 


-4.8 


1.50 


±0.47 


+3.8 
-3.6 


-9.0 




±0.6 


±1. 


.6 


+0.13 
-0.03 


+4.2 
-9.8 



Table 4. Measured £^ enslty and systematic uncertainty breakdown for the dijet data. The system- 
atic uncertainties marked with a * are uncorrelated between \rj\ bins. 



{ Jl d J ){UE) Stat - E i E 2 M x M 2 M 3 * P J Total 



r-^){MB) 



[7oj 


0.0 


- 0.8 


2.95 


±0.64 


0.8 


- 1.6 


2.81 


±0.57 


1.6 


- 2.4 


2.61 


±0.55 


2.4 


- 3.2 


2.43 


±0.52 


3.2 


-4.0 


2.21 


±0.53 


4.0 


-4.8 


2.00 


±0.49 



70 70 70 70 70 70 70 70 



+1.1 +1.5 +0.27 in A +1.6 +4.3 

-1.3 -1.6 -0.28 ^t-J-^ _i.3 _4.2 

+1.7 +0.64 +0.10 +0.25 , A n +1.3 +5.1 

-1.6 -0.69 -0.11 -0.26 _i.3 _5.i 

+1.2 +0.08 +0.43 in r +0.98 +3.9 

-1.2 -0.08 -0.47 _o.92 -3.9 

+1.4 +2.3 _i_o a +0.80 +3.8 

-2.1 -2.5 _o.37 -4.2 

+1.7 _ 9 9 +1.5 in a +0.46 +4.1 

-1.6 z - z -1.6 =c<J- i ± _o.l7 -4.6 

+1.1 _ 9 a +0.19 in n +0.13 +4.1 

-1.0 -0.20 =c<J-y -0.03 -4.7 



Table 5. Ratio of measured E T cnsl y for the dijet data to that for the the minimum bias data, and 
systematic uncertainty breakdown. The systematic uncertainties marked with a * are uncorrelated 
between |fy| bins. 
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1 dN cvt 


Stat. 


E? E 2 


Mi 


M 2 


MS 


Pi 


Total 


JV C vt dSE T 


[GeV] 


[GeV- 1 ] 


[%] 


[%] [%] 


[%] 


[%] 


[%} 


[%] 


[%] 


0-2 


0.161 


0.35 


4-9 o 
-2.1 


=F 2.1 




+ 0.38 


±0.07 


I o n 
-\~0.\J 

-3.1 


2-4 


0.0835 


0.32 


-r±.o 
-1.4 


T 1-2 




+ 0.21 


=F4.5 


J-A R 

-4.8 


4-6 


0.0533 


0.40 


-0.39 


+ 0.90 




+ 0.16 


+3.8 


-3.9 


6-8 


0.039 


0.46 


— U.ZU 

+0.00 


± 0.23 




± 0.04 


+0.44 


in ao 
-t"U.Do 

-0.71 


8-12 


0.0271 


0.49 


—0.57 
+0^56 


± 1.4 




± 0.24 


±2.2 


4-9 7 

-2.7 


10 1 R 
1Z — ID 


U.Ul 1 1 


n K7 

U.D 1 


-1.3 
+ 1.5 


± l.o 




4- n 


4-Q ^ 
±o.O 


+4.3 
-4.2 


16 - 20 


0.0121 


0.67 


-2.5 
+2.5 


± 2.4 




± 0.43 


±3.6 


+5.1 
-5.1 


20 - 30 


0.00619 


0.75 


-4.6 
+4.8 


± 4.1 




± 0.73 


±4.2 


+7.7 
-7.5 


30 - 40 


0.00226 


1.2 


-7.5 
+8.4 


± 7.9 




± 1.4 


±5.5 


+13 
-12 


40 - 50 


0.000855 


1.9 


-10 
+ 12 


± 8.7 




± 1.5 


±8.1 


+17 
-16 


50 - 60 


0.000321 


2.5 


-13 
+ 15 


± 13 




± 2.3 


±10 


+22 
-21 



Table 6. Measured j^- ggjg; and systematic uncertainty breakdown for the minimum bias data 
in the region 0.0 < \t]\ < 0.8. 



[GeV] 
0-2 

2-4 

4-6 

6-8 

8-12 

12 - 16 

16 - 20 

20 - 30 

30 - 40 

40 - 50 

50 - 60 



1 dN cvt 
N cvt dT,E T 

[GeV" 1 ] 



Stat. 

[%1 



0.123 


0.38 


0.092 


0.30 


0.0588 


0.35 


0.0425 


0.42 


0.0296 


0.44 


0.0198 


0.51 


0.0137 


0.58 


0.00726 


0.67 


0.00268 


1.1 


0.000951 


1.7 


0.000339 


2.4 



E\ 

[%] 

+4.7 

-4.8 

+2.5 

-2.9 

+0.55 

-0.85 

-0.29 
+0.21 

-1.0 
+1.1 
-2.2 
+2.0 
-3.9 
+3.8 
-7.3 
+7.6 
-13 
+14 
-18 
+21 
-22 
+30 



E2 



Mi 


M 2 


M 3 


Pi 


[%] 


[%] 


[%} 


[%] 


+ 3.0 


+ 0.46 


T 1.1 


±0.64 


+ 1.6 


+ 0.25 


+ 0.62 


+5.3 


+ 0.65 


T 0.1 


+ 0.25 


+9 


± 0.07 


± 0.01 


± 0.03 


+3.8 


± 0.76 


± 0.12 


± 0.29 


±1.3 


± 1.9 


± 0.29 


± 0.73 


±5.7 


± 2.1 


± 0.32 


± 0.80 


±8.3 


± 4.7 


± 0.72 


± 1.8 


±10 


± 7.3 


± 1.1 


± 2.8 


±9 


± 13 


± 1.9 


± 4.9 


±10 


± 17 


± 2.6 


± 6.5 


±15 



Total 

[%] 

+5.7 

-5.8 

+6.1 
-6.2 

+9.0 
-9.1 
+3.8 
-3.8 
+1.9 
-1.9 
+6.4 
-6.4 
+9.4 
-9.4 
+14 
-13 
+19 
-17 
+27 
-25 
+38 
-32 



Table 7. Measured ]v^T3&F anc ^ systematic uncertainty breakdown for the minimum bias data 
in the region 0.8 < \r)\ < 1.6. 
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1 dN cvt 
N avt dT,E T 


Stat. 


E\ 


E12 


Mi 


M 


2 


Mf 


Pi 


Total 


[GeV] 


[GeV- 1 ] 


[%] 


[%] 


[%] 


[%] 


[% 


»] 


[%] 


[%] 


[%] 


0-2 


0.0980 


0.41 


-LA A 

-4.3 


~~ 


~~ 


=F 


1.2 


T 6.8 


±1.1 


-1-8 9 

-8.2 


2-4 


0.0931 


0.32 


-1-9 ^ 

i-z.o 
-2.6 






=F 


0.48 


T 2.7 


=F2.0 


J_A 9 

-4.3 


4-6 


0.0639 


0.36 


-0.76 


_ 


_ 


=F 


0.18 


T 1.0 


T6.2 


Ifi Q 

-6.4 


6-8 


0.0460 


0.45 


—0.15 
+0.05 






± 


0.27 


± 1.6 


=F4.8 


-ru.i 
-5.1 


8-12 


0.0323 


0.46 


—0.61 

+0^55 






± 


0.26 


± 1.5 


=F0.75 


-1-1 K 

-|--L.O 

-1.8 


10 1 fi 
1Z — ID 


n noi p. 

U.UZ10 


n ^/i 

U.04 


-1.6 
+1.5 






± 


u.oo 


4- 1 Q 

± i.y 


4-0 1 
±Z.l 


+3.2 
-3.3 


16 - 20 


0.0149 


0.62 


-2.9 
+2.7 






± 


0.59 


± 3.4 


±3.5 


+5.6 
-5.7 


20 - 30 


0.00792 


0.68 


-5.6 
+5.6 






± 


1.2 


± 7.1 


±6.8 


+ 11 
-11 


30 - 40 


0.00290 


1.1 


-11 
+11 






± 


2.6 


± 15 


±12 


+22 
-22 


40 - 50 


0.000977 


1.8 


-15 
+17 






± 


3.2 


± 18 


±14 


+29 
-28 


50 - 60 


0.000312 


2.5 


-19 
+24 






± 


4.3 


± 25 


±14 


+37 
-34 



Table 8. Measured j^- ggjg and systematic uncertainty breakdown for the minimum bias data 
in the region 1.6 < \t]\ < 2.4. 



[GeV] 


1 dN cvt 


Stat. 

[%] 


Ef 
[%} 


E2 
[%} 


Mi 
[%} 


M 2 
[%] 


M 3 
[%} 


Pi 

[%] 


N evt dZE T 

[GeV" 1 ] 


0-2 


0.0840 


0.48 


+6.9 
-3.1 








T 5.3 


±0.72 


2-4 


0.0920 


0.33 


+3.1 
-1.6 








T 1.7 


Tl.5 


4-6 


0.0660 


0.46 


+0.65 
-0.41 








± 0.46 


T5.3 


6-8 
8-12 


0.0485 
0.0343 


0.47 
0.50 


-0.16 
+0.07 
-0.46 
+0.25 








± 0.47 
± 0.33 


^3.5 
T0.57 


12 - 16 


0.0232 


0.57 


-1.5 
+0.6 








± 1.1 


±3.1 


16 - 20 


0.0160 


0.61 


-3.7 
+1.6 








± 2.7 


±3.3 


20 - 30 
30-40 


0.00829 
0.00290 


0.72 
1.2 


-7.8 
+3.6 
-14 

+7 








± 4.6 
± 7.6 


±5.3 
±6.9 


40 - 50 


0.000908 


1.6 


-22 
+12 








± 12 


±11 


50 - 60 


0.000281 


3.2 


-27 
+14 








± 15 


±7.3 



Total 

[%] 

+8.8 

-6.2 

+3.9 

-2.8 

+5.3 

-5.3 

+3.5 

-3.5 

+0.86 

-0.95 

+3.4 

-3.7 

+4.6 

-5.7 

+8 

-11 

+ 12 

-17 

+20 

-27 

+22 

-32 



Table 9. Measured ^ 

iv cv t ai^hji 

in the region 2.4 < \rj\ < 3.2. 



and systematic uncertainty breakdown for the minimum bias data 
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1 dNcvt Q-f n + pe 
AT evt dT,E T SlaZ - ^1 


£ 2 


Mi M 2 M| 




Total 


[GeV] 


[GeV" 1 ] [%] [%] 


[%] 


[%] [%] [%] 


[%] 


[%] 



n 

u - 


Z 


n ns^7 

U.UOO / 


n ap. 


2 - 


4 


0.0931 


0.33 


4 - 


6 


0.0720 


0.38 


6 - 


8 


0.0530 


0.44 


8 - 


12 


0.0367 


0.47 


12 


- 16 


0.0237 


0.54 


16 


- 20 


0.0154 


0.66 


20 


- 30 


0.00704 


0.77 


30 


- 40 


0.00183 


1.4 


40 


- 50 


0.000385 


2.4 



-H 7.3 T4.0 Tl-5 tH 

Iju 4.9 T 2.1 ±0.64 1^ 

±1:1 1-6 ±0.18 T 4.1 itjj 

-oil 0-33 ±0.57 T4.5 1^ 

+oil "0-98 ±1-6 +H 

+ii "3-3 ± 1-1 ±2.3 +3.2 

_6.4 - - ±2.6 ±5.3 1™ 

7_i;jj -13 ±5.7 ±9.5 I** 

;H -25 ±7.5 ±13 t\\ 

~\\ -33 ± 20 ±13 +H 



Table 10. Measured dZET an< ^ systematic uncertainty breakdown for the minimum bias data 



in the region 3.2 < |?7| < 4.0. 



Y,E T f3tat ^ M\ ^ 2 M 3 Pl Total 

[GeV] [GeV- 1 ! [%] [%] [%] [%] [%] [%] [%] [%] 

9.0 — — =F 1-7 T0.37 ±1% 

11° 51 ± - 58 =Fi-i ±1:1 

18:?? 1-5 ±0.54 T5.3 1^ 

+°- 06 0.16 T 0.14 T4.5 +4 - 5 



- 


2 


0.0818 


0.54 


2 - 


4 


0.102 


0.42 


4 - 


6 


0.0777 


0.44 


6 - 


8 


0.0577 


0.48 


8 - 


12 


0.0397 


0.49 


12 


- 16 


0.0239 


0.58 


16 


- 20 


0.0135 


0.75 


20 


- 30 


0.00464 


0.96 


30 


-40 


0.000642 


2.3 



-0.04 u - iu T -r^- u -4.5 

-1.5 ± 0.50 T0.26 t°i 



-6.0 ± 0.33 ±3.6 tj A 5 

~H -12 - - ±1.6 ±9.1 11° 
~ 8 - 7 -22 ± 2.5 ±14 +17 



+8.9 ^ - 1 - -28 

+11 -39 ±3.4 ±24 



Table 11. Measured ]v^3&F an< ^ systematic uncertainty breakdown for the minimum bias data 



in the region 4.0 < \rj\ < 4.8. 
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1 dN cvt 
7V C vt dT,E T 


Stat. 


^2 


Mi 


M 2 


MS 


P2 


J 


Total 


[GeV] 


[GeV" 1 ] 


[%] 


[%] [%] 


[%] 


[%] 


[%] 


[%] 


[%} 


[%] 


0-2 


0.0998 


2.1 


-4.4 


T 2.1 




T 0.37 


T5.1 


— 1.2 
+L2 


4-7 d 

-7.4 


2-4 


0.0870 


1.5 


-4-3 4 
-3.7 


=F 1.4 




T 0.24 


±2.6 


—0.61 
+0.35 


-1-4 8 
-5.0 


A 

4i — 


U.U 1 01 


1 A 
1.44 


+2.0 
-1.8 


+ U.U 1 




-+• U.1Z 


-1-9 7 


-0.38 
+0.08 


+3.7 
-3.6 


6-8 


0.0598 


1.6 


+0.46 
-0.59 


± 0.05 




± 0.01 


Tl-8 


-0.34 
+0.14 


+2.5 
-2.5 


8-12 


0.0409 


1.7 


-1.9 
+1.9 


± 1.1 




± 0.20 


T2.0 


-0.23 
+0.25 


+3.4 
-3.4 


12 - 16 


0.0224 


2.2 


-4.9 
+5.3 


± 2.5 




± 0.45 


T0.15 


+0.14 
+0.21 


+6.3 
-6.0 


16 - 20 


0.0121 


3.0 


-7.7 
+7.0 


± 4.0 




± 0.70 


±6.0 


+1.0 
-0.8 


+ 11 
-11 


20 - 30 


0.00428 


4.3 


-11 
+13 


± 6.5 




± 1.1 


±3.9 


+5.4 
-4.5 


+16 

-15 


30 - 40 


0.000939 


8.7 


-16 
+18 


± 10 




± 1.8 


±16 


+13 
-11 


+31 

-28 



Table 12. Measured and systematic uncertainty breakdown for the dijet data in the 

region 0.0 < |r?| < 0.8. 





1 dN cvt 


Stat. 


E\ E2 


Mi 


M 2 


M| 


P2 


J 


Total 


N cvt dZE T 


[GeV] 


[GeV- 1 ] 


[%] 


[%} [%} 


[%] 


[%] 


[%] 


[%} 


[%] 


[%] 


0-2 


0.0680 


2.2 


+9.7 
-9.6 


T 4.0 


T 0.62 


T 1.5 


±7.7 


-0.79 
+0.89 


+13 

-13 


2-4 


0.0844 


1.5 


+7.1 
-7.4 


T 2.8 


T 0.44 


T 1.1 


±4.8 


-0.80 
+0.81 


+9.2 
-9.4 


4-6 


0.0804 


1.4 


+4.1 
-4.4 


T 1.7 


T 0.25 


T 0.64 


±3.7 


-0.83 
+0.91 


+6.0 
-6.2 


6-8 


0.0695 


1.4 


+0.9 
-1.6 


T 0.48 


T 0.07 


T 0.19 


±2.0 


-0.62 
+0.69 


+2.7 
-3.0 


8-12 


0.0476 


1.6 


-3.4 
+2.7 


± 1.3 


± 0.20 


± 0.49 


T0.86 


-0.01 
+0.09 


+3.5 
-4.1 


12 - 16 


0.0252 


2.0 


-8.9 
+8.9 


± 3.6 


± 0.56 


± 1.4 


+"4.2 


+0.69 
-1.1 


+ 11 
-11 


16 - 20 


0.0129 


2.8 


-13 
+15 


± 6.0 


± 0.92 


± 2.3 


^4.8 


+ 1.3 
-2.3 


+ 17 
-16 


20 - 30 


0.00429 


4.2 


-19 
+23 


± 10 


± 1.6 


± 3.9 


±3.5 


+4.2 
-4.2 


+26 
-23 


30 - 40 


0.000785 


9.5 


-26 
+32 


± 16 


± 2.5 


± 6.1 


±6.9 


+ 13 
-10 


+40 
-35 



Table 13. Measured -^—^ and systematic uncertainty breakdown for the dijet data in the 

region 0.8 < | ^7 1 < 1.6. 
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1 dN cvt 
AT . Ay 7?^ 


Stat. 


El 


E 2 


Mi 


M 2 


Ml 


P2 


J 


Total 


[GeV] 


[GeV- 1 ] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%} 


[%] 


[%] 


[%] 


0-2 


0.0604 


2.4 


+7.6 
-7.5 


— 


— 


T 1-6 


T 9.2 


=F8.2 


+0.01 
+0.52 


+15 

-15 


2-4 


0.0859 


1.6 


+6.1 
-6.0 


— 


— 


=F 1-2 


T 6.6 


±2.6 


-0.43 
+0.51 


+9.6 
-9.5 


4-6 


0.0857 


1.5 


+3.4 
-3.3 






T 0.70 


T 4.0 


±2.1 


—0.47 
+0.27 


+5.8 
-5.8 


6-8 


0.0715 


1.6 


+0.8 
-1.3 






T 0.24 


T 1.3 


±2.5 


-0.40 
+0.22 


+3.3 
-3.5 


8-12 


0.0485 


1.7 


-2.7 
+2.2 






± 0.46 


± 2.6 


±0.51 


-0.63 
+0.33 


+3.9 
-4.2 


12 - 16 


0.0256 


2.3 


-7.3 
+7.4 






± 1.4 


± 7.9 


±3 


-0.30 
+0.10 


+ 12 
-11 


16 - 20 


0.0124 


3.3 


-12 
+12 






± 2.3 


± 13 


±1.2 


+0.75 
-0.51 


+ 18 
-18 


20 - 30 


0.00385 


4.8 


-17 
+18 






± 4.0 


± 22 


±1.3 


+4.7 
-3.9 


+30 

-29 


30 - 40 


0.000641 


11 


-20 
+25 






± 6.3 


± 36 


+-2.0 


+ 17 
-14 


+48 
-45 



Table 14. Measured and systematic uncertainty breakdown for the dijet data in the 

region 1.6 < < 2.4. 



YjE^ 
[GeV] 


1 dN cv t 


Stat. 

[%] 


Ef 
[%} 


E 2 

[%] 


Mi 
[%] 


M 2 
[%] 


M| 
[%] 


P2 

[%} 


JVcvt dSE T 

[GeV" 1 ] 


0-2 


0.0565 


2.9 


+n 

-5 








T 2.6 


+^6.0 


2-4 


0.0892 


1.7 


+8.5 
-3.8 








T 1-8 


±0.53 


4-6 


0.0905 


1.7 


+5.1 
-2.5 








T 0.9 


±2.7 


6-8 


0.0748 


1.7 


+1.1 
-1.0 








T 0.03 


±1.3 


8-12 


0.0495 


1.8 


-4.6 
+1.1 








± 1.3 


Tl.2 


12 - 16 


0.0253 


2.6 


-12 

+6 








± 3.0 


±1.5 


16 - 20 


0.0111 


3.5 


-18 
+10 








± 4.8 


±2.1 


20 - 30 


0.00298 


5.8 


-25 
+15 








± 7.8 


±1.2 



Total 

[%] 

+13 

-9 

+8.9 

-4.5 

+6.1 

-4.2 

+2.4 

-2.3 

+2.8 

-5.3 

+7 

-13 

+12 

-19 

+18 

-27 



Table 15. Measured dZET anc ^ systematic uncertainty breakdown for the dijet data in the 



region 2.4 < |?7| < 3.2. 
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1 iWovt 
AT JV U 1 

iVevt aZj^T 


Stat. 


E\ 


E2 


Mi 


M 2 


M| 


Pi 


Total 


[GeV] 


[GeV- 1 ] 


[%] 


[%} 


[%} 


[%] 


[%] 


[%] 


[%] 


[%] 


0-2 


0.0784 


2.5 


+8.3 
-7.3 


12 






T 5.1 


T-6.6 


+17 
-11 


2-4 


0.105 


1.7 


+6.7 
-6.2 


9.5 






T 2.9 


±1.7 


+12 
-7.3 


4-6 


0.0999 


1.6 


+2.8 
-3.1 


3.9 






T 0.8 


±2.0 


+5.5 
-4.1 


6-8 


0.0756 


1.7 


-1.4 
+0.3 


-2.0 






± 1.3 


±0.47 


+2.2 
-3.3 


8-12 


0.0433 


2.0 


-7.1 
+5.6 


-10 






± 4.5 


±0.46 


+8 
-13 


12 - 16 


0.0172 


2.9 


-15 
+ 15 


-21 






± 8.8 


±0.25 


+17 
-27 


16 - 20 


0.00601 


4.7 


-19 
+25 


-27 






± 13 


T0.76 


+29 
-36 



Table 16. Measured dEEh anc ^ systematic uncertainty breakdown for the dijet data in the 



region 3.2 < |?7| < 4.0. 



YjE'y 


1 oWcvt 


Stat. 


4 


E2 


Mi 


M 2 


M l 


P2 


Total 


7V cv t cfE&r 


[GeV] 


[GeV- 1 ] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%} 


[%} 


[%] 


0-2 


0.0915 


2.5 


+6.2 
-5.4 


16 






=F 0.63 


^6.2 


+ 18 

-9 


2-4 


0.139 


1.6 


+3.4 
-3.5 


8.5 






T 0.12 


±1.3 


+9.3 
-4.1 


4-6 


0.113 


1.6 


-0.13 
-0.51 


-0.33 






± 0.38 


±0.38 


+ 1.7 
-1.7 


6-8 


0.0726 


1.8 


-3.6 
+2.7 


-9.0 






± 0.88 


T0.18 


+3.4 
-9.8 


8-12 


0.0313 


2.5 


-7.9 
+8.4 


-20 






± 1.6 


±4.7 


+ 10 

-22 


12 - 16 


0.00775 


4.4 


-11 
+13 


-28 






± 2.6 


±2.2 


+ 14 

-31 



Table 17. Measured dSET anc ^ systematic uncertainty breakdown for the dijet data in the 
region 4.0 < |?7| < 4.8. 
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